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In this note we discuss metastability in a long-but-finite range disordered model for the glass 
transition. We show that relaxation is dominated by configuration belonging to metastable states 
and associate an in principle computable free-energy barrier to the equilibrium relaxation time. 
Adam-Gibbs-like relaxation times appear naturally in this approach. 

I. INTRODUCTION 

Metastable states and configurational entropy, are central concepts in our understanding of low temperature relax- 
ation of disordered and glassy systems [1] . In presence of time scale separation between "fast" and "slow" degrees of 
freedom relaxation is often described as a walk between metastable states, intended as regions of configuration space 
where the fast degrees of freedom reach a state of quasi equilibrium before the relaxation of slow degrees of freedom 
can effectively take place. According to the Adam and Gibbs theory [2] the configurational entropy i.e. the logarithm 
of the multiplicity of metastable states is related to the relaxation time by an inverse proportionality relation. 

The relation between metastability and relaxation is well understood in infinite range spin glasses with "random 
first order transition" [3] . The dynamics of these systems is exactly described by the equations of Ideal Mode Coupling 
Theory [4]. The Mode Coupling transition, spurious in liquid systems, is here driven by the existence of metastable 
states capable to trap the system for time scales diverging with the systems size. The divergence is of course a Mean 
Field artifact, but features of mean- field metastable states could well have a reflex in finite dimensional glassy systems 
[5]. 

A well known problem in the study of metastable states is that while they are easy to define within mean field theory 
their definition becomes elusive in finite dimensions. This problem has been studied in large detail in the context 
of systems with first order phase transition, where the problem of metastability, can be considered as satisfactorily 
understood [6,7]. Unfortunately this is not the case in glassy systems, where despite the appeal of a description of 
dynamics in terms of metastable states, only very rare contributions have tried to clarify the meaning of metastability 
beyond the phenomenological level [8]. 

The "random first order transition scenario" has led to phenomenological attempts to treat the problem of glassy 
relaxation, and ergodicity recovery in analogy with kinetics of ordinary first order phase transition [9-11]. Liquids at 
low temperature would appear as "mosaic states" , in which locally frozen variables' configurations could be composed 
combinatorially on a large scale. The typical size of the rigid regions could be computed in a nucleation theory, 
with the configurational entropy playing the role of an ergodicity restoration bulk driving force, competing with some 
postulated interface tension. These developments stimulated first principle calculations in microscopic disordered 
models. In ref. [12] and later in [13] a statistical description of the ergodic state and a computation of the glassy 
coherent length below T c was proposed through the asymptotic analysis of a Landau-like glassy effective free-energy 
derived from microscopic models. Unfortunately in that papers it was not possible to make an explicit connection 
between the computed free-energy barrier and the relaxation time of the system. Scope of this letter is to discuss 
this connection. In order to do that, we start from an analysis of glassy relaxation based on separation of time scales. 
We argue that finite dimensional relaxation is dominated by metastable states that can be characterized along the 
classical lines of Lebowitz and Penrose (LP) [14], first proposed to describe metastable phases of matter in presence 
of first order phase transitions. According to LP, metastable states can be considered as constrained equilibrium 
ensembles with: 1) an order parameter is homogeneous on a suitable mesoscopic length scale 2) a large time life and 
3) a very small probability of return once abandoned. In analogy with the work of LP, we use in our analysis models 
with long-but-finite range Kac kind of interactions, which in our case have a disordered character. These offer the 
possibility of studying finite dimensional effects in an expansion around mean field, and the local mean-field character 
of correlations, postulated in the phenomenological mosaic description, appears as a consequence of the range of 
interaction [15]. 
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II. THE MODEL 



Let us consider a spherical Kac p-spin glass model [12] defined for real spins Oi on the d-dimensional hypercubic 
lattice A of linear size L and Hamiltonian 

H(a,J) = - Jii-i P <rii---<ri P (!) 

where the couplings Ji 1 ..-i p are i.i.d. Gaussian variables with zero average and variance 

E{Jl.. ip ) = \l pd £ V(7Ki - k\) ■ ■ ■ ^\i p k\) (2) 
fceA 

where p is an integer p > 3 and -0 ( ) , is a non-negative integrable function verifying the normalization J d d x ip(\x\) = 
1. With this choice, the couplings J^.-.i are sensibly different from zero only if all pairs of variables \ii — i TO |^ 7 _1 
V7 7^ m = 1, so that only variables that are at distances \i— j] < -f^ 1 effectively interact. The effective interaction 
range -f' 1 will be assumed to be large throughout the paper. We partition the lattice in boxes B x of a coarse graining 
length size I and impose a local spherical constraint J2 ieB a 1 = ^ ^ or au x - We are interested to the regime where 
the three defining lengths are supposed, as in LP, to verify the relation 

logL«£« 7 - 1 «L (3) 

and for dcfiniteness we will have in mind the situation where I — j~ s for some S £ (0, 1). The model is chosen in 
such a way to reduce to the usual mean-field spherical p-spin model in the regime 7L sa 1 . Let us recall the physics 
of this case [16] that will be useful in the finite 7 case. On lowering the temperature from the paramagnetic region 
one encounters two transitions. There is a first transition at the Mode Coupling temperature T c where ergodicity 
is broken. Below T c an exponential number of ergodic components Af « e ^s(T) dojninateg the thermodynamics. 
The configurational entropy S(T) decreases as a function of the temperature, until it becomes zero at a "Kauzmann 
temperature" Tk < T c , where a second transition is met. 
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FIG. 1. Qualitative behavior of the mean field potential. From top to bottom the T = T c , Tk < T < T c and T = Tk- The 
difference of free-energy between the two minima is the product of the configurational entropy ad the temperature. 

This structure of states can be studied through the so-called "quenched glassy potential" that considers the free- 
energy of a system a constrained to have a fixed global overlap q — ^ i a i<3* with a typical reference equilibrium 
configuration a* [17]. This function has the shape depicted in figure 1, in the region Tk < T < T c one finds a two 
minimum structure analogous to the free-energy as a function of the order parameter for mean-field systems with a 
first-order phase transition. The value of the overlap in the low minimum and high minima represent respectively 
the typical overlap among random configurations in different and the same ergodic components. Both minima are 
associated to the exponential multiplicity of ergodic components, and their difference gives directly the configurational 
entropy times the temperature <&(<7ea) - $(0) = TE(T). 
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III. DYNAMICS AT FINITE 7. 



We now address the question of restoration of ergodicity for finite 7 and Tk < T < T c . The relaxation of model 
(1) can be conveniently studied through Langevin dynamics 



doi , . dH 



(4) 



where fx x (t) are Lagrange multipliers that enforce in average the local spherical constraint at each time and the rji 
are white noise variables with correlations E(r)i{t)r)j{s)) — 2T5ij5(t — s). An insight on the dynamics in the regime 
(3) can be obtained studying the Kac limit 7^0 after the thermodynamic limit. Through standard dynamical 
techniques -Martin-Siggia-Rose formalism [18] or dynamical cavity method [19]-, it is possible to conclude [20] that 
the finite time local equilibrium spin-spin correlation function C^\t) = jj J2 ieB 0»(i)0i(O) on the scale £, tends to 
a space homogeneous function, C(t), that for temperatures larger then the Kauzmann temperature Tk verifies the 
usual equilibrium mean field equation for the spherical p-spin model [16] 



dsC p -\t-s) 



dC(s) 
ds 



(5) 



This equation predicts ergodicity breaking below T c . The equilibrium correlation function in this regime approaches 
exponentially fast on a scale tmf a n on zero limit lim^oo C(t) = q^A- Conversely, for any finite 7 > 0, for finite or 
infinite L, the system should be able to escape from the metastable states and lim^oo lim^^oo C x L,1 \t) coincide with 
the equilibrium correlation function jj ^2 ieB (&i) 2 = 0. In the regime (3), the behavior of C x j) (t) can be schematically 
depicted as in figure 2. As a function of time, the small 7 correlation function remains close to the mean field curve 
up to some characteristic time diverging when 7 _1 — > 00 and it decays to zero on longer scales. Since for small 
enough 7 this characteristic time is arbitrary larger than the mean field relaxation time tmf, one finds that the system 
equilibrates in the region of phase space where the local overlaps with the initial condition are close to q&A before 
relaxing further. Moreover, one can expect that on time scales of the order of tmf, the system has homogeneous 
correlations on the scale of the box size £, with typical fluctuation from cell to cell of (i) — C(t) scaling as l/l d / 2 . In 
turn, assuming independence of the fluctuations beyond some finite length, using extreme statistics one can estimate 

the maximum deviation max x [C2(t) — C(t)] to be at most of order This, according to our hypothesis, can be 

made arbitrarily small. The previously discussed properties are enough to characterize the regions explored on scales 
smaller then as metastable states in the Lebowitz-Penrose sense. 
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FIG. 2. Qualitative behavior of the correlation function as a function of time for different values of 7. 

The role of the "order parameter" , homogeneous in space, is played by the overlap of the configuration of the system 
with the initial condition <r(0). Formally, for any typical equilibrium configuration er*, we can define metastable 
states, as equilibrium ensembles where the configurations a are such that in each box B x the overlap with the initial 
configuration q x (a,a*) = -p X^es^ a i a i 1S suitably close to q B A- More precisely, one defines an interval I = [<j_, 1], 
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with g_ < qEA and the set 1Z = {a\q x (a, a*) E I}. The metastable state is the restricted ensemble 

Mrc(o-) = -^-exp(-pH(a))l n (a). 

As in the LP case, the definition of metastable state makes sense for a whole range of values of g_, such that a local 
fluctuation of the overlap larger then g_ is with high probability pushed back to qEA- However, in order for the 
definition to be useful, the value of q_ should be fixed a posteriori in such a way that a fluctuation larger then that 
size should be amplified, and eventually run the system out of the metastable state with finite probability. We can 
now use this analysis to estimate the relaxation rate of our system as a decay rate of the metastable state. We can 
then proceeds as in LP, and argue that for the local Langevin dynamics we are using, the relaxation rate A^ = -^j 
is proportional to A^ 7 ^ « ; |F L where dIZ is the set of configuration a belonging to 1Z, such that q Xo (a, a*) = q- for 
at least one point xq. A^ 7 ' is in principle a random variable that depends on the initial configuration a* and on 
the quenched random couplings J^,...,^. However the relaxation rate should be self-averaging with respect to both 
sources of noise, and its typical value can be computed considering the average of its log 

log A^ = £(log Z dn - log Z n ) (6) 

where E stands for the canonical average over a* and the quenched average over the distribution of the disorder 
J. Eq. (6) is the main result of this paper. It connects the relaxation rate to a free-energy barrier computable in 
principle in purely static terms. 

IV. RELATION WITH THE EFFECTIVE POTENTIAL. 

For fixed disorder and given finite volume we can define an "effective potential" as the free-energy cost to keep the 
overlap with an equilibrium configuration a* in each of the boxes B x fixed to specified values q x : 

eM-^W a ,,j[q x }) = £ J] O - Qx)e~ pH ^ (7) 

^ ax 

where in the expression, we have anticipated the scaling of the free-energy with the interaction range 7 _1 [21]. Is is 
clear how Zg-ji and Z-jz can be related to W^*^^], in fact, for S = 1Z, dIZ, with transparent notation 

Z S = j ' <Dq x e-$ W "' Aq * ] . (8) 

By the previous discussion, the dominant profile in Z-r is estimated as the constant q x = qEA for all x. For this 
profile, the extensive part of W a * t j[q x — qEA] is a self-averaging quantity: V^ba] « E(W a * t j[q x = qEA])- 

For generic overlap profiles q x , we can consider V[q x ] — E(W a * ,j[q x ]), and ask if the knowledge of this quantity is 
enough to reconstruct the relaxation rate for small 7. 

The restricted partition function Zqh is: 

Zen = Y.I Vq x e-^'^\ 

x "'feo=9-; q^el Va;} 

Each term of the sum on xq will be dominated by a single profile, minimizing Wo-*^^] in presence of the constraints. 
These will be profiles with a localized excitation around xq , corresponding to a finite free-energy cost with respect to 
the uniform state q x — qEA- 

For fixed a* and J, different values of x could give rise to different free-energy cost and different profile shapes. 
However, we could argue again that on the mesoscopic scale 7 _1 almost all sites £0 would give the contributions of the 
same order and these dominate the relaxation rate. This is a stronger hypothesis than the other ones underlying our 
reasoning. It implies that the incipient spatial heterogeneities on a scale 7" 1 responsible for relaxation, have equal 
probability of appearing in any point of space, independently of the reference configuration a* and the realization of 
the quenched disorder J. Without that, the computation of A^ should start from an evaluation of (6). This is in 
principle possible and when done could allow the validation or refutation of the hypothesis. 



5 



Previous work [10,12,13] was implicitly based on that strong hypotheses. Within that framework, if we denote by 
q* the profile maximizing V[q x ] for xq = 0, the relaxation rate will be given by 



and V[q*] will provide a Landau free-energy functional that we can use in a nucleation-like theory. It is clear that 
restricted to profiles constant in space, q x = q, for small 7 we have: V(q x = q) = W a *.j\q x = q] w (L-f) d Q(q), 
where $ is the mean-field potential. In order to compute the free-energy barrier within this self-averaging setting, 
one should consider spherically symmetric solutions with boundary conditions limi^i^oo q x = qEA- Formulated in this 
way the problem reduces to the computation of a critical droplet as suggested by phenomenological theories [9,11]. 
According to conventional nucleation theory [6], close to Tk where the two minima are close to degeneracy, one can 
use the "thin wall approximation" and evaluate the interface tension through the ID instantonic solutions connecting 
the two minima. This competes with a bulk driving term given by the free-energy difference between the minima, 
which in our case is the mean-field configurational entropy. This first principle formulation confirms the prediction of 
phenomenological theories, that the configurational entropy provides the bulk driving force to ergodicity restoration. 
The actual computation of V for fixed profile can be performed through a modified version of the replica method [17], 
needed to average over the spin reference configuration a* and the quenched disorder J. This has been done for the 
model (1,2) in ref. [12]. Without entering in the details of the discussion of that paper we just mention that within 
a specific replica ansatz it is possible to find a finite surface tension is found for T approaching Tk- This implies a 
glassy coherent length behaving as R g ~ S(T) -1 and an Adam-Gibbs free energy barrier AV ~ ■ This result 

is supported by a more elaborated ansatz in [13]. As discussed in [12], and in view of possible non self-averaging 
effects, these solution could in fact just provide approximations to the true barrier. The "right" replica solutions could 
change the exponent d — 1 appearing in the relation. However, the general Adam-Gibbs structure of the formula, 
that only depends on the fact that the bulk ergodicity restoring force is the configurational entropy, should remain 
unchanged. 



Summarizing, we have presented an argument that identifies as Lebowitz-Penrose metastable states the dominant 
configurations of relaxing glassy systems below the Mode Coupling transition. This allows to compute the relaxation 
time in purely static term through the evaluation of the probability of equilibrium fluctuations. The relation between 
relaxation time and free-energy barrier is summarized in eq. (6) which is the main result of this letter. We have 
also discussed how previous results in the literature relate implicitly to a self-averaging hypothesis. If this holds, the 
relevant fluctuations have the shape of overlap droplets, and the bulk driving force for ergodicity restoration is provided 
by the configurational entropy. In agreement with phenomenological theories the relation of the relaxation time with 
the configurational entropy follows a generalized Adam-Gibbs equation. The argument has been presented in the case 
of disordered spin model, where playing with the interaction range one can have time scales as much separated as 
wanted and a mesoscopic scale for which correlations are homogeneous on short time scales. On the other hand our 
argument is actually only based on the hypothesis that time separation is such to allow partial equilibration of the 
fast degrees of freedom for fixed values of the slow one, and as such it should be generically valid when this condition 
is respected. 
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